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In this paper we use the Nuclear Magnetic Resonance (NMR) to write eletronic states of a ferro- 
magnetic system into a high-temperature paramagnetic nuclear spins. Through the control of phase 
and duration of radiofrequency pulses we set the NMR density matrix populations, and apply the 
technique of quantum state tomography to experimentally obtain the matrix elements of the system, 
from which we calculate the temperature dependence of magnetization for different magnetic fields. 
The effects of the variation of temperature and magnetic field over the populations can be mapped in 
the angles of spins rotations, carried out by the RF pulses. The experimental results are compared 
to the Brillouin functions of ferromagnetic ordered systems in the mean field approximation for two 
cases: the mean field is given by (i) B — Bo + AM and (ii) B = Bo + AM + A'M 3 , where Bo is the 
external magnetic field, and A, A' are mean field parameters. The first case exhibits second order 
transition, whereas the second case has first order transition with temperature hysteresis. The NMR 
simulations are in good agreement with the magnetic predictions. 
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INTRODUCTION 

Magnetism and magnetic materials are among the 
main branches of research in modern Condensed Mat- 
ter Physics. Magnetic materials have been extensively 
studied for decades by a variety of experimental tech- 
niques, including Nuclear Magnetic Resonance (NMR) 
In a standard NMR experiment, a sequence of ra- 
diofrequency (RF) pulses is applied to a sample, which 
may or may not be subject to a static external magnetic 
field (for the last case, zero- field NMR). Typical infor- 
mation obtained in such NMR experiments are: local 
magnetic fields and moments, spin and charge distribu- 
tions, magnetic anisotropy, relaxation dynamics, etc. For 
reviews of NMR on magnetic materials see Ref. [l[ and 
references therein. 

For the last ten years NMR has been established also 
as an useful tool for Quantum Information Processing 
(QIP), and several quantum algorithms and protocols 
have been tested in the period. The reason for this is the 
fact that RF pulses are equivalent to unitary transforma- 
tions, which are needed for QIP, and can be controlled 
with great precision in NMR experiments, allowing the 
manipulation of the quantum states of a system (HQ and 
generation of protocols to process quantum information 
li-SI . Particularly interesting is the use of NMR-QIP to 
simulate quantum systems [7H11[. 



In an usual quantum simulation performed by NMR, 
the quantum dynamics of a given system is emulated by 



mapping the Hamiltonian of the system into the NMR 
Hamiltonian [1, [l2-15|. In this work, we use NMR to 
write electronic states of a ferromagnetic ordered mate- 
rial into the paramagnetic nuclei density matrix. This 
is made by first calculating the trace distance between 
the respective density matrices for various temperatures, 
and then working out the pulse sequences to achieve the 
correct eletronic state. The results show how well nuclear 
states can be used as a quantum memory and, although 
the present study does not exploit the influence over the 
eletronic magnetization curve, it is our belief that such 
studies could bring to NMR a new useful way to study 
the magnetization dynamics of materials. Some possible 
applications arc discussed on the conclusions. 



MEAN-FIELD MODEL OF MAGNETIC 
SYSTEMS 



Magnetic systems are usually modeled using the 
Heisenberg Hamiltonian which, in the mean-field approx- 
imation, can be simplified to [l|: 

H = -gfi B B-8, (1) 

and the electronic magnetization is given by (for the z- 
componcnt of spins): 

M = gfi B tr[p(T)S z ] = g fx B 8 B s {x) , (2) 

where p(T) stands for the thermal equilibrium density 
matrix, x = gps B S/ksT and Bg(x) is the Brillouin 



2 



function. In a paramagnetic isolated spin system, the 
field B is just the external field B$. In the mean field 
approximation, on the other hand, exchange interaction 
is parameterized by an effective field B = B +XM, where 
A = 2(5 — l) 2 zJ ex /g 2 p 2 B is called 'mean-field parameter', 
g is the Lande g— factor, z the number of first neighbors, 
p b the Bohr magneton, J ex the exchange-parameter and 
Eq.([2]) must be solved self-consistently [I]. One of the 
main aspects of the mean-field approximation is that it 
predicts magnetic ordering below a critical temperature 
T c = g 2 p 2 B S (S + 1 ) A/3 ks- From this relation we can 
observe that the critical temperature gives us a measure 
of the exchange parameter A. 

In particular, the Brillouin function for a spin S = 3/2 
is given by: 



B 3 /2(a0 = ^{2coth(2x)-icoth(|)| 



(3) 



Equations ([2]) and ([3]) will be used to compare the mean- 
field prediction for an ordered magnetic system with 
the NMR results of 23 Na in a room temperature liquid- 
crystal sample. 



NMR TWO QUBIT SYSTEMS 

A nuclear system composed by two interacting nuclei, 
A and B, with 1 = 1/2 spins in a static magnetic field 
is the prototype of a two-qubit NMR quantum computer 
0, The Hamiltonian of such a system is: 



/2.1/2 



La 1 z,A 



1t^l b Iz,b + 27rJ I z a I z 



(4) 



where ul a{b) are the respective Larmor frequencies, and 
2ttJ is the coupling constant. Such a Hamiltonian de- 
scribes a four-level system, for instance, coupled H and 
13 C spins in a molecule of Chloroform, under a static 
magnetic field. Alternatively, a four-level system can be 
described by a NMR system composed by a spin 3/2 nu- 
cleus in a static magnetic field and in the presence of a 
local electric field gradient. In such case, the Hamiltonian 
is given by [l|: 

Hujq 



U 



3/2 



-hw L i z 
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(5) 



where ujq is the quadrupole coupling constant. This sys- 
tem can be used to emulate any two-qubit quantum sys- 
In NMR systems the Zeeman energy 



tern m 



levels are very small if compared to the thermal energy 
at room temperature. This means that at room tempera- 
ture the NMR density matrix can be written as [El. [2#]26| : 



p=-I + eAp, 



(6) 



where I is the 2" x 2™ identity matrix, n is the number 
of qubits, Ap = I z is the deviation density matrix and 
e = hu L /2 n k B T~ lCT 5 . 



Unfortunately, such state given in Eq.© is inadequate 
for quantum computing proposes, since for that it is nec- 
essary a pure inicial state [5|. However, the NMR ability 
for manipulating spins states resulted in a method for 
creating states isomophic to a pure state, named pseudo- 
pure states [27], H3- Such pseudo-pure states (PPS), typ- 
ically have the form: 



P 



(1-e) 



epi 



(7) 



There are some different methods to create these states 
||. In this work we have used the time-averaging 
method; the basic pulse sequences for generating pseu- 
dopure states in a / = 3/2 quadrupole system are given 
in Refs. [Il-[31| • 



Upon the application of a sequence of radio-frequency 
pulses representing an unitary transformation, U, Eq.© 
transform accordingly: 



UpU^ 



(1-e) 



I + eUp^ 



(8) 



From the application of a sequence of such pulses, by 
setting the pulses duration, phase and amplitude, a very 
fine control over the density matrix population and co- 
herences can be achieved, and it is possible to generate all 
two-qubit computational base states, also superposition 
and entangled states 22, 2t| . It is important to note that 
these entangled states are called pseudo-entangled states 
because e ~ 10~ 5 , which makes Prm S always separable 
even when pi is an entagled state [32| . 



WRITTING FERROMAGNETIC STATES INTO 
NMR STATES 

Equation |J2J) can be used to write ferromagnetic elec- 
tronic states into a high-temperature paramagnetic nu- 
clear spin system. As a prototype we consider a coupled 
two-spin 1=1/2 system. The initial state p\ corresponds 
to the pseudo-pure state [13,[2§| |00}(00|. The sates are 
labeled as 1 00) , |01), |10), and |11), for increasing order of 
energy, in accordance to the current literature of NMR 
quantum information, where each spin 1/2 represents a 
qubit. 

Considering rotation angles Of and Of , of the two spins 
(A and B) over the directions x and y by the operator 

Rx,y(0 x ,0f ,0y ,6y) = R X {8 X ) R X (0 X ) Ry(9 y ) Ry(0y), 

where the rotation matrices are given by: 



R i 



(gi I B cos ( 
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Rj(6f) =I A ®I B cos^ 



> If sin 



(9) 



(10) 
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and j = x, y, we arrive at the following populations: 
Poo = i(l + cos6£ cos 0£)(1 + cos 6»f cosf?^), 



00 



PlO 



Pll 



i = |(1 + cos 0^ cos0;f)(l-cos0f cos0f) 



i(l-cos6# cos0^)(l + cos0f cos0^), 



(11) 



cos0^ cos0^) (1 



-cos0f cos0f ) . 



^ and 0^, one can set the 



Therefore, by controlling 
NMR levels populations. 

In order to obtain the energy levels populations from 
the density matrix p(T) of Eq.® we use the concept of 
trace distance [331 ] : 



D 



tr(|p-p'|) 



(12) 



where p represents the density matrix of the NMR sys- 
tem and p' the target density matrix. For numerical pur- 
poses, for each temperature we seek for a pulse sequence 
for which D < 10 -3 . This establishes a relationship be- 
tween the rotation angles and the temperature T . With 
this kind of mapping we are able to manipulate the pop- 
ulation of the pseudo-pure states in order to mimic the 
corresponding electronic population at a given tempera- 
ture. Note that, since we are dealing with pseudo-pure 
states, their population can be manipulated in such a way 
that any temperature can be mimetized including T = 
K. 

The above description using two spin 1/2 is very conve- 
nient because it provides an analytical description for the 
rotation operators and populations. However, the imple- 
mentation can also be achieved in the spin 3/2 quadrupo- 
lar system described by the Hamiltonian in Eq.([5]). In 
this case, the states 1 00) , |01), 1 10) , and |11) can be asso- 
ciated to the four energy levels of the spin 3/2 system and 
pseudo-pure states and rotation operators that act inde- 
pendently in each qubit can be built, in complete analogy 
with the spin 1/2 system. A convenient manner of creat- 
ing the pseudo-pure states and the qubit selective rota- 
tion in this system is using a numerically optimized pulse 
sequence, named Strongly Modulated Pulses (SMP)[34|- 
371 ] . Therefore, because in our experiments we used the 
quadrupolar spin 3/2, all rotation operators were imple- 
mented with the SMP technique, and the experimental 
density matrices were reconstructed with the quantum 
state tomography method described in Ref. [13] ■ 

Our experiments have been carried out on 23 Na (/ = 
3/2) nuclei dissolved in a lyotropic liquid crystal, de- 
scribed by the Hamiltonian in Eq.©. The sample was 
prepared with 20.9 wt% of sodium dodecyl sulfate (95% 
of purity), 3.7 wt% of decanol, and 75.4 wt% deuterium 
oxide, following the procedure described elsewhere [38j|. 
The 23 Na NMR experiments were performed using a 
9.4 T - VARIAN INOVA spectrometer using a 7-mm 
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FIG. 1: Theoretical results (continuous lines) and NMR im- 
plementation (open symbols) for a simple ferromagnet at 
Bo = (open squares) and Bo — 6 T (open circles). 



solid-state NMR probe head. The quadrupolc coupling 
luq /27r was about 16 kHz. 



Second order phase transition 

Ferromagnetic systems described by Eq.([2|), where 
B = Bo + AM, g = 2 and z = 6, have a second or- 
der phase transition, e.g., the solution of these mean- 
field model agrees very well with the solution for the 
Ising spin chain with infinite size. From Bo = it is 
possible to obtain the critical temperature of the sys- 
tem T c = g 2 p? B S(S + 1) A/3 k B . The different tempera- 
tures were implemented by manipulating the populations 
of the energy levels, and such manipulations were done 
using the radio-frequency (RF) pulses. By minimizing 
the trace distance, Eq. (|12[) . between the density matrix 
elements in Eq. lfTTj) and the ones of Eq.([2]), we could 
map the temperatures to be emulated into rotation an- 
gles of the RF pulses that originates differences among 
the populations. This made possible the emulation of 
the temperature. For the electronic system considered 
system here, the coherences do not exist. Thus, to emu- 
lated the temperature correctly in the NMR system, we 
used in the SMP implementation with temporal averag- 
ing techniques [B[ in order to cancel this elements of the 
density matrix. 

In FigfT]the theoretical calculation and the NMR ex- 
perimental results for the magnetization as a function of 
temperature curve of a simple ferromagnet are presented. 
The continuous line represents the Brillouin function for 
Bo = and B Q = 6 T. The open symbols are the NMR 
results, in good agreement with the ferromagnet curve 
prediction. 
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First order phase transition 



CONCLUSIONS 



Interesting additional features appear in mean field ap- 
proximation if we add an extra M 3 term in the effective 
magnetic field, B = B + AM + A'M 3 . This mean-field 
model presents a first order phase transition together 
with a thermal hysteresis, as observed in some magnetic 
systems 3^, 4(|. This hysteresis appears due to a differ- 



ence on the behavior of the Gibbs free energy between 
heating up (superheating) and cooling down (supercool- 
ing) the system (40l | . It is known from Landau theory that 
both A and A' parameters rule the first order transition; 
especially the critical temperature T c and the superheat- 
ing and supercooling temperatures (those that limit the 
thermal hysteresis). In other words, the occurence of the 
phase transition and the degree of the hysteresis depend 
on A and A'. Thus, being A = 2(g — I) 2 zJ ex / g 2 n 2 B with 
g = 2, z = 6, and J ex = 0.5 x 10 _19 fce, we choose A' <C A 
(typically A'/A ~ ICT 2 ) to obtain T c = 83 K when cool- 
ing down and T c = 60 K when heating up, which opti- 
mize the view of the phase transition and the thermal 
hysteresis. 

In Fig|2]we show the NMR implementation of such sys- 
tem. The continuous lines represent the Brillouin func- 
tions and the open circles the NMR implementation. For 
Bq = 0, the first order phase transition appears. Also in 
Fig |2] it can be seen the heating up of the system and 
the cooling down. Our NMR implementation described 
quite well the thermal hysteresis expected theoretically. 
For Bq > B c , it is theoretically expected that the jump 
on the magnetization disappears, and our NMR imple- 
mentation also described it quite well. 
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FIG. 2: Theoretical results (continuous lines) and NMR im- 
plementation (open symbols). For Bo — 0, the thermal hys- 
teresis appears and is well described in the NMR system. For 
Bo > B c (open circles) , the jump on the magnetization disap- 
pears and it is also well written on the NMR system. The open 
triangles are for temperature sweeping up and open squares 
are for temperature sweeping down. 



In this work we successfully described the temperature 
and magnetic field effects over two distinct ferromag- 
netic systems through NMR quantum information pro- 
cessing techniques: the ferromagnetic ordering described 
by two mean-field models. Differently from the NMR 
works found in literature, we directly emulated the den- 
sity matrix of the ferromagnetic system by establishing a 
relationship between NMR spin rotations and tempera- 
ture. This was done by minimizing the trace distance be- 
tween the NMR and the target system density matrices. 
The NMR experiments correctly exhibit first and second 
order phase transitions, as well as thermal hysteresis. As 
a proof of principle, the model systems used here have 
analytical solution, so the main point of this article is 
to show that high temperature nuclear spin systems can 
be manipulated to behave as ordered electronic spin sys- 
tems. Given that, we believe that this technique can be 
extended to simulate other interesting magnetic phenom- 
ena, for instance magneto-caloric effect, study of critical 
exponents, quantum phase transitions, etc. Particularly 
interesting would be the study of the environment effects 
over the magnetization, by looking at the magnetization 
phenomena. 
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